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ABSTRACT 

We show that images of TeV blazars in the GeV energy band should contain, along with point-like 
sources, degree-scale jet-like extensions. These GeV extensions are the result of electromagnetic cas- 
cades initiated by TeV 7-rays interacting with extragalactic background light and the deflection of the 
cascade electrons/positrons in extragalactic magnetic fields (EGMF). Using Monte-Carlo simulations, 
we study the spectral and timing properties of the degree-scale extensions in simulated GeV band 
images of TeV blazars. We show that the brightness profile of such degree-scale extensions can be 
used to infer the lightcurve of the primary TeV 7-ray source over the past 10 7 yr, i.e. over a time scale 
comparable to the life-time of the parent active galactic nucleus. This implies that the degree-scale 
jet-like GeV emission could be detected not only near known active TeV blazars, but also from "TeV 
blazar remnants" , whose central engines were switched off up to ten million years ago. Since the 
brightness profile of the GeV "jets" depends on the strength and the structure of the EGMF, their 
observation provides additionally information about the EGMF. 

Subject headings: gamma rays: galaxies - galaxies: active - galaxies: jets - methods: numerical 



Introduction. Significant progress in understanding the 
activity of blazars, i.e. active galaxies with relativistic 
jets aligned with the line of sight, was achieved with the 
start of operation of the Fermi telescope. The combi- 
nation of data from Fermi in the 0.1 - 10 GeV energy 
band and from ground based 7-ray telescopes like HESS, 
MAGIC and VERITAS in the 100 GeV - 10 TeV band 
provides detailed simultaneous spectral and timing infor- 
mation for t he most extreme r epresentatives of the blazar 
population (jAbdo et alJl2009t l. 

The TeV 7-ray flux from distant blazars is signif- 
icantly attenuated on the way from the source to 
the Earth by pair production on the i nfrared/optical 
extra g alactic backgr o und l ig ht (EBL) ( Knciske ~et alj 
20041: iStecker et al l 120061 : iFranceschini et all 12008 
Primack et al. 1 120081) . TeV 7-rays that are absorbed 
on the way from the primary 7-ray source initiate 
electromagnetic cascades in the intergalactic space. 
Charged component of the electromagnetic cascade is 
deflected by the EGMF. Potentially observable effects 
of such electromagnetic cascades in the EGMF i nclude 
the delayed "echoes" of m ulti-TeV 7-ray flares (jPlagal 
1995s IMurase et all |2008) and the appearance of ex- 
tended emission around initially point-like 7-ray sources 
(lAharonian et al.lll99l iNeronov fc Semikod 120071 120091: 
Dola g et all I200a lElviv et al.ll2009l) . 

TeV 7-ray emission from blazars is believed to be rela- 
tivistically beamed into a narrow cone (jet) with an open- 
ing angle 8 jet ~ r _1 - 5° [r/10], where V is the bulk 
Lorentz factor of the 7-ray emitting plasma. Blazars are 
a special type of 7-ray emitting AGN for which the angle 
between the line of sig ht (LOS) and the jet axis, o bs, is 
0obs < e je t, see Fig. Q] (lUrrv et al.l ll99l. 

In general, the number of blazars with a given jet-LOS 



.©, 



jel 



1 £ 


/©ex, 


■^L___]6obs 







Fig. 1. — Geometry of the propagation of direct and cascade 
7-rays from the source (on the left) to the observer (on the right). 

misalignment angle is expected to scale as dN/d9 h s ~ 
6* bs in the range < 9 ^ s < Oj C t- This means that most 
of the observed TeV blazars should have 9 b s ~ @jet, 
rather than # b s — 0. Consequently, the TeV 7-ray 
emission pattern is not symmetric with respect to the 
axis source-observer. In sources with o b s ~ ©jet, most 
multi-TeV 7-rays are emitted preferentially on one side 
of the LOS, as it is shown in Fig. U As a result, the 
extended emission from the cascade initiated by the ab- 
sorption of TeV 7-rays in interactions with EBL photons 
should appear as one-sided extension next to the primary 
point source, rather than as an extende d axially symmet- 
ric halo discussed in previous studies (lAharonian et al] 
1991IN eronov fc Sem ikoz 2007, 2009; IDolae et al.H 2009: 
Elviv et al.ll2009D . 

In what follows we discuss the spectral and timing 
properties of such one-sided cascade extensions in the 
0.1-1 GeV energy band images of TeV blazars. Our 
study is based on two independent Monte Carlo codes 
for the modeling of 7-ray induced electromagnetic cas- 
ca des in the in t ergala ctic space, which were introduced 
bv lDolag et all pul) and lElviv et~aH ((2009T) . 

Basic formulae. Before presenting our numerical re- 
sults, we discuss the basic physics of the phenomenon in 
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a simplified picture. In particular, we assume that the 
electromagnetic cascade consists only of two steps and 
replace probability distributions by their means. Then 
the mean free path of VHE 7-rays through the EBL can 
be approximated by 
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where the numerical factor k — k(E Jo ,z) ~ O(l) in- 
cludes the uncertainties of the EBL modeling. Pair pro- 
duction on EBL reduces the flux of 7-rays from the source 
by a factor exp [— r(i? 7o )], where r ~ D/D 7 is the opti- 
cal depth with respect to pair production and D is the 
distance to the source. 

Electron-positron pairs created in interactions of multi- 
TeV 7-rays with EBL photons produce secondary 7- 
rays via inverse Compton (IC) scattering on Cosmic Mi- 
crowave Background (CMB) photons. Typical energies 
of the IC photons reaching the Earth are 
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where ecmb = 6 x 1CP 4 eV is the typical energy of CMB 
photons. 

Deflections of e + e~ pairs produced by the 7-rays, 
which were initially emitted away from the observer, can 
redirect secondary photons toward the observer. This ef- 
fect leads to the appearance of exten ded emission around 
an initially point source of 7-rays ( Neronov fc Semikozl 
[2007t IDolag et al.ll2009l ; lElviv et al.|[2009D . 

In the absence of perfect alignment of the jet axis with 
the line of sight, the extended cascade emission might 
be strongly asymmetric. It might appear as a jet-like 
feature next to the primary 7-ray source. The maximal 
angular size of this jet-like feature can be estimated as 
the size of the projected 7-ray mean free path as 
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if Dy < D (i.e. r > 1). If r(E la ) < 1, the cascade 
emission from the TeV 7-ray beam can extend to very 
large angles 6 cx t .max 

— 7T/2. 

The jet-like extended emission can be observed only if 
deflections of the cascade e + e~ pairs are sufficiently large 
to redirect the cascade emission toward the observer. If 
the correlation length of the EGMF is larger than the 
electron cooling distance 
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where <jt denotes the Thomson cross section and E/cmb 
the energy density of the CMB p hotons, then the de- 
flecti on angle can be estimated as (jNeronov fe Semikozl 
2009) 
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with Rl as the Larmor radius of electrons and positrons. 

If the correlation length Xb of the EGMF is much 
smaller than the electron cooling distance D e , the de- 



flection angle can be estimated using the diffusion ap- 
proximation as 
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If the EGMF is weak, electron/positron trajectories 
are not strongly deflected during one cooling time and 
thus secondary cascade 7-rays are emitted within a cone 
with opening angle of order 0(S). In this case only a part 
of the cascade emission could be observed. If the mean 
free path of the primary 7-rays is much shorter than the 
distance to the source, the angular extension could be 
estimated from the simple geometrical consideration of 
Fig. Q]as 

sin(e oxt (B)) = ^-sinS. (7) 

Otherwise, the angular size of the source is found from 
the sum of the angles of triangle with vertices at the 
source, at the pair production point and at the position 
of the observer (see Fig. [IJ as 



e ext (s) = 5- 
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Most of the known TeV blazars have moderate dis- 
tances, so that t(E J0 = 1 TeV) < 1. In this case, 
6 C xt,max ~ tt/2 and Qext(B) = 5 - #obs- A measure- 
ment of O cx t <§C S xt,max thus provides a measurement 
of 6 and, in this way, gives a constraint on the parameters 
of the EGMF, i.e. B and Xb- 

The difference in the path length between the direct 
and cascade 7-rays leads to a significant time delay of 
the cascade emission signal. For a given jet misalignment 
angle o bs, the time delay of emission coming from the 
direction 9 away from the source is 
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Comparing this time scale with the typical time scale of 
AGN activity, Tagn ~ 10 7 Y T : one sees that degree-scale 
extended emission in the GeV energy range depends on 
the TeV 7-ray luminosity of the blazar integrated over 
its life time. 

Results of numerical modeling. To model the asym- 
metric extended emission from the 7-ray initiated elec- 
tromagnetic cascade in intergalactic space, we have 
extended the two Mo nte-Carlo codes in troduced by 
IDolag et all (f2009T) and lElviv et al.l (|2009h : The codes 
follow now the three-dimensional trajectories of individ- 
ual cascade particles moving through the EGMF. The 
turbulent component of t he EGMF has been calcula ted 
following the algorithm of IGiacalone fc Jokipiil ((1994). 

To produce an image of the 7-ray induced electro- 
magnetic cascade, as it would be detected by a 7-ray 
telesco pe, we use the algorithm described by lElviv et al.l 
(2009). We have verified that the results obtained using 
the two different codes are compatible with each other. 
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Fig. 2. — E > 1 GeV band images of the sky region around TeV 
blazars with jets inclined at o b s = 0°, 9 b s = 3°, 8 b s = 6° and 
Sobs = 9° (left to right). The jet opening angle is 0j ct = 3°, the 
EGMF strength is B = 10~ 16 G. The spacing of the coordinate grid 
is 2°, the color scale is logarithmic in surface brightness: Yellow 
corresponds to the maximal surface brightness, black corresponds 
to the surface brightness less than 10 — 3 of the maximal value. 

We record positions and directions of all secondary 7- 
rays which cross a sphere of the radius R = D around the 
source. We choose the directions of primary 7-rays to be 
distributed within a cone with an opening angle Oj C t- We 
consider a primary 7-ray beam with a Gaussian profile, 
so that the probability for a primary 7-ray to have a 
direction misaligned by an angle with respect to the 
jet axis is p(6) ~ exp ( — 9 2 /0j Ot ). 

For simplicity, we consider a monochromatic primary 
7-ray beam with all the primary 7-rays having the same 
energy E J0 — 1 TeV. This is sufficient to demonstrate the 
existence of the effect discussed here for the first time, 
namely degree-scale jet-like extensions in Fermi images 
of TeV blazar s. The EBL backgrou nd is taken from the 
calculations of lKneiske et al.l (|2004f) . We fix the distance 
to the source as D — 400 Mpc, so that r(E l0 ) ~ 1. The 
EGMF is chosen to have a correlation length of the order 
of several Mpc, with its power spectrum sharply peaked 
at the wavenumber k ~ I Mpc . Our results could 
be generalized in a straightforward way to the case of 
an arbitrary primary 7-ray spectrum, arbitrary distance 
to the source and different EBL models, when consider- 
ing extended emission from particular TeV blazars with 
known TeV band spectra and known redshift. 

Figure [5] shows the effect of the misalignment of the 
primary 7-ray beam with the LOS on the morphology 
of the extended emission. The left panel of the Figure 
corresponds to the situation o h s = 0, which is equivalent 
to the axia lly-symmetr i c case considered by iDolag et al.1 
(2009) and lElviv et al.l (|2009D . An axially-symmetric ex- 
tended " halo" around the primary point source is clearly 
visible. The other panels of the Figure show the cases 
of a jet with opening angle Oj 0t = 3° misaligned by the 
angles 9 b s = 3°, 6° and 9°, respectively. It is clear that 
the misalignment of the jet axis with the line of sight 
leads to the appearance of an extended jet-like feature 
on one side of the source. The ratio of the point source 
flux to the flux of the extension grows with the increase 
of the misalignment angle o h s - 

The angular extension of the cascade emission depends 
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Fig. 3. — E > 1 GeV band images of the sky region around TeV 
blazars with ©j et = 9 Q b s = 3° for different values of the EGMF 
strength. From left to right: 10" 17 G, lCT 16 G, 10" 15 G, lO" 14 G. 
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Fig. 4. — E > 1 GeV band images of the sky region around 
a TeV blazar with 0j ct = S bs = 3° at different times following 
instantaneous injection of 1 TeV 7-rays at the source. From left to 
right: images in time intervals < idol < 10 5 yr, 10 5 yr< tjgj < 
10 6 yr, 10° yr< i del < 3 X 10 6 yr and 3 X 10 6 yr< t del < 10° yr 
after the outburst. 

on the strength of the EGMF as long as the trajectories of 
e + e~ pairs are not completely randomized. The morpho- 
logical properties of the jet-like emission are practically 
independent from the properties of the EGMF, when the 
EGMF strength is such that the deflection angle 6 > 2ir, 
cf. Eqs. (0) and ([6]). Figure [3] shows the growth of the 
source extension with the increase of the EGMF strength. 
For magnetic fields stronger than B ~ I0~ 15 G, the size 
of the extended source reaches ten(s) of degrees. In this 
case, the extended source could significantly contribute 
to the diffuse 7-ray background. 

Cascade emission coming from regions with angular 
distance 8 > 1° to the primary source is delayed by 
2deiay ~ 10 5 -10 7 yr compared to the direct emission from 
the source (see Eq. (flQ|) ). This means that "echos" from 
periods of enhanced activity of the source (e.g. an en- 
hanced accretion rate following major merger episodes), 
which happened all along the life time of an AGN some 
time T ago, could enhance the flux at the distance 
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9 ~ 1.7° [T/10 6 yr] [(9 ohs + 9 je t)/5°] from the source. 

Figure 14] shows a time sequence of E > 1 GeV band 
images of the sky region around a TeV source at differ- 
ent times after a short episode of TeV 7-ray emission. 
One can clearly see that the emission at large angular 
distances is delayed by up to 10 7 yr. 

The flux coming from the region at an angular distance 
9 from the point source is proportional to the source flux 
averaged over the period T<j e i a y Therefore it is possi- 
ble that GeV 7-rays are detectable today from an AGN 
which was active some 10 7 yr ago, but is at present not 
active anymore. In this case a GeV source would be clas- 
sified as "unidentified": The parent AGN (a) could not 
be identified as an AGN in the optical, X-ray and TeV 
7-ray bands or (b ) the GeV source is displaced from the 
position of the parent AGN. The characteristic feature 
of such an unidentified "AGN remnant" is the absence 
of counterparts at lower energies. If the GeV 7-rays are 
produced by e + e~ pairs deposited in the intergalactic 
medium by primary TeV 7-rays, the only energy loss 
mechanism for the pairs is inverse Compton scattering 
on CMB photons. 

Discussion. The presence of extended jet-like emission 
at degree scales should be a generic feature of GeV band 
images of TeV blazars. The total flux of the jet-like ex- 
tended source is proportional to the source luminosity in 
the TeV energy band. Taking into account the fact that 
TeV blazars have hard 7-ray spectra, the primary source 
luminosity in the TeV band could be much larger than 
its GeV luminosity, so that the overall extended source 
luminosity could be higher than the primary source lu- 
minosity in the GeV band. This means that the best 
candidates for the search of extended emission are TeV 
blazars with hard intrinsic spectra. 

This does not automatically mean that the extended 
emission should be readily detectable in Fermi images 
of TeV blazars. In spite of the larger luminosity, the 
extended source flux might be suppressed if the EGMF 
is strong enough to randomize the trajectories of e + e~ 
pairs before they loose their energy to the GeV band 
via inverse Compton emission. The maximal possible 
suppression of the extended source flux is by a factor 
9r t 2 ~ 100. 

Another potential problem for the detection of jet-like 
extended emission next to TeV blazars is that the ex- 
tended source has to be identified on top of the diffuse 
7-ray background (which is composed of the Galactic and 
extragalactic diffuse emission). The minimal detectable 
flux for extended sources increases roughly as 9 1 ' 2 , where 
9 is the angular length of the jet-like extended source. In 
this respect, sources at larger distances, for which the 
jet-like extensions appear more compact, are better can- 
didates for the search of extended emission in the Fermi 
energy band. 
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Finally, the detectability of extended emission close to 
TeV blazars strongly depends on the angular resolution 
of the LAT telescope. At low energies, E~ ~ 0.1 GeV, the 
LAT angular resolution is relatively poor, 6*psf — 10°. It 
is clear that only very large angular size jet-like exten- 
sions with an angular diameter 9 ~ #psf could be de- 
tected. However, the detectability of such large extended 
sources would be complicated by the high level of dif- 
fuse sky background within the ~ 10° region around the 
source. At the same time, the size of the PSF decreases 
to 9psf < 1° above GeV energies. This dramatically im- 
proves the sensitivity of the telescope for the search of 
extended emission: Extensions of much smaller angular 
size could be detected on top of a strongly reduced back- 
ground. This favors the search of extended emission at 
energies above ~ 1 GeV. 

The detection of degree-scale jet-like extensions in 
Fermi images of TeV blazars would provide new insight 
into the physics of 7-ray emission from blazars. It might 
provide information about the history of activity of the 
supermassive black hole in the blazar central engine over 
a time span comparable to typical AGN life times. The 
study of the morphology of the extended emission would 
make a measurement of the opening angles 0j C t of blazar 
jets possible. It would also give information about the 
strength and correlation length of the EGMF. 

To summarize, we have shown that GeV band images 
of TeV blazars should possess degree-scale jet-like ex- 
tended features. These features trace the direction of 
the TeV 7-ray beam emitted by the blazar. They are 
produced as results of electromagnetic cascades initiated 
by TeV 7-rays interacting with EBL photons. We have 
performed Monte Carlo simulations of three-dimensional 
electromagnetic cascades developing in the EGMF. Us- 
ing these Monte Carlo simulations, we have derived the 
properties of the GeV jet-like extended emission near 
TeV blazars. We have investigated the dependence of 
the characteristics of the jet-like extended sources (the 
angular size, the brightness profile) on the strength of 
the extragalactic magnetic fields and on the opening an- 
gle and orientation of the primary TeV 7-ray beam from 
the blazar. We have also demonstrated that the 7-ray 
signal in the jet-like extended emission is delayed up to 
10 7 years compared to the direct 7-ray signal from the 
primary point source. The very long time delay of the 
cascade emission means that the extended GeV source 
could be detected next to a blazar which is no longer 
active as blazar. 
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